Feedback Cooling of the Normal Modes of a Massive Electromechanical System to 

Submillikelvin Temperature 
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We apply a feedback cooling technique to simultaneously cool the three electromechanical normal 
modes of the ton-scale resonant-bar gravitational wave detector AURIGA. The measuring system is based 
on a dc Superconducting Quantum Interference Device (SQUID) amplifier, and the feedback cooling is 
applied electronically to the input circuit of the SQUID. Starting from a bath temperature of 4.2 K, we 
achieve a minimum temperature of 0.17 mK for the coolest normal mode. The same technique, imple- 
mented in a dedicated experiment at subkelvin bath temperature and with a quantum limited SQUID, could 
allow to approach the quantum ground state of a kilogram-scale mechanical resonator. 

PACS numbers: 42.50.Lc, 04.80.Nn, 05.40.-a, 45.80.-r 



The prospect of observing quantum behaviour and in- 
vestigating decoherence in macroscopic mechanical res- 
onators 01 0] requires cooling to ultralow temperatures, 
such that the thermal energy becomes comparable to the 
quantum energy. In addition to conventional refrigera- 
tion, this goal requires the use of feedback cooling tech- 
niques [3D, unless the frequency is as large as several hun- 
dred MHz. Various authors have recently reported ad- 
vances in cooling mechanical resonators, using either op- 
tomechanical JiLlI IS II IE IIE 111 or electromechan- 
ical [T2l techniques, and implementing either active 
feedback n, pi |q, ]31 or dynamical back-action effects 
IE & fiol tlu lrlll . These experiments involved mostly 
nanomechanical or micromechanical resonators, and have 
demonstrated cooling capability down to mK temperatures 
|@]. Cooling of a gram-scale optical spring resonator to 
a few mK has been demonstrated as well using the tech- 
niques developed for interferometric gravitational wave 
(GW) detectors 0]. In this Letter, we show that very ef- 
ficient cooling can be achieved even in much larger sys- 
tems, exploiting the techniques developed for resonant- 
mass GW detectors, based on electromechanical transduc- 
ers coupled to Superconducting Quantum Interference De- 
vice (SQUID) sensors. In this case, one can take advan- 
tage of the larger quality factor achievable in macroscopic 
systems with respect to micromechanical ones. As exper- 
imental demonstration, we simultaneously cool the three 
electromechanical normal modes of the ton-scale resonant- 
bar GW detector AURIGA. Our cooling technique is based 
on an electronic feedback directly applied to the input cir- 
cuit of the SQUID. 

AURIGA represents the state-of-art in the class of res- 
onant GW detectors [13], and has been in continuous op- 
eration from year 2004, searching for galactic astrophysi- 
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FIG. 1: (Color online) Scheme of the gravitational wave detector 
AURIGA. The system comprises three coupled resonators with 
nearly equal resonant frequency of about 900 Hz: the first longi- 
tudinal mode of the cylindrical bar, the first flexural mode of the 
mushroom shaped resonator, which is also one of the plates of 
the electrostatic capacitive transducer, and the low-loss electrical 
LC circuit. The electrical current of the LC resonator is detected 
by a low noise dc SQUID amplifier. 



cal events in collaboration with a world network of detec- 
tors 03]. It is located in Padua (Italy) and is based on a 
2.2 x 10 3 kg, 3 m long bar made of a low loss aluminium 
alloy (A15056), cooled to liquid helium temperatures. The 
fundamental longitudinal mode of the bar, sensitive to 
gravitational waves, has an effective mass M=l.l x 10 3 
kg and a resonance frequency w B /27r = 900 Hz. Accord- 
ing to the equipartition theorem, the rms amplitude of the 
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where k B is the Boltzmann constant and T is the temper- 
ature. For the AURIGA bar, the rms thermal motion is 
x rms = 4 x 10~ 17 m at T = 4.2 K. The bar resonator 
motion is detected by a displacement sensor with a sen- 
sitivity of order several 10~ 20 m/\/Hz over a ~ 100 Hz 
bandwidth. This sensitivity is accomplished by a multi- 
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mode resonant capacitive transducer [15] combined with 
a very low noise dc SQUID amplifier flql (Fig. [TJ. In 
this scheme, the bar resonator is coupled to the funda- 
mental flexural mode of a mushroom-shaped lighter res- 
onator, with 6 kg effective mass and the same resonance 
frequency. As the mechanical energy is transferred from 
the bar to the lighter resonator, the motion is magnified by 
a factor of roughly 15. A capacitive transducer, biased with 
a static electric field of 10 7 V/m, converts the differential 
motion between bar and mushroom resonator into an elec- 
trical current, which is finally detected by a low noise dc 
SQUID amplifier through a low-loss high-ratio supercon- 
ducting transformer. The transducer efficiency is further 
increased by placing the resonance frequency of the elec- 
trical LC circuit close to the mechanical resonance frequen- 
cies HI, at 930 Hz. 

The detector can then be simply modelled as a system 
of three coupled resonators: its dynamics is described by 
three normal modes at separate frequencies, each one being 
a superposition of the bar and transducer mechanical res- 
onators and the LC electrical resonator lfl3ll . As seen from 
the SQUID sensor, each mode k (k= 1,2,3) is modelled as 
a RLC series electrical mode with an effective inductance 
L k , capacitance C k and resistance R k (Fig. 12). The to- 
tal inductance L k includes also the input inductance L in 
of the SQUID amplifier. Around the resonance frequency 
of each mode uj k /2ir, with io k = (L k C k ) , the com- 
plex impedance of the circuit is expressed by Z k (co) = 
R k + iivL k + 1/ (iujC k ). We point out that, although the 
modes appear as purely electrical as seen from the SQUID, 
their dynamics actually includes the full motion of both 
mechanical resonators. In fact, the current in the electri- 
cal circuit is linearly related with the mechanical motion 
of bar and transducer, and the effective impedance param- 
eters L k , C k , R k of each normal mode are determined by 
the mechanical and electrical parameters of all resonators 
lfl3ll . Thus, cooling the normal modes of the system im- 
plies cooling the motion of the mechanical resonators. The 
effective impedance parameters L k , C k , R k are experimen- 
tally estimated by measuring the current I s = V ca \/Z k 
in response to a calibration voltage signal V ca \, injected 
through a small calibration coil in series to the SQUID 
input coil (see Fig. 0*. The inductances L k of the three 
modes are respectively 1.66 x 10~ 4 H, 1.23 x 10~ 5 H and 
8.12 x 10~ 6 H. The mode resonance frequencies LO k /2n 
are 865 Hz, 914 Hz, 953 Hz, and the quality factors Q k = 
uj k L k /R k are 1.2 x 10 6 , 0.88 x 10 6 , and 0.77 x 10 6 . 

The effective temperature T k of each mode is propor- 
tional to the mean square current (I k ) induced by thermal 
fluctuations driving the mode, according to the equiparti- 
tion theorem: 

(It) = kj f±. (1) 

To reduce the mode temperature, we implement an elec- 
tronic feedback cooling (Fig. O, or cold damping, tech- 
nique [ 17]. The SQUID output voltage V = AI S is passed 
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FIG. 2: (Color online) The feedback cooling scheme. The k-th 
normal mode is approximated, around its resonance frequency, 
by a series-RLC circuit. In this representation, different modes 
should be thought of as being in parallel with each other. The 
dc SQUID is represented as current amplifier. The electronic 
feedback cooling is obtained by sending back a current Id phase 
shifted of 7r/2 with respect to I s - 



through a passive low-pass filter D (uj) with cut-off fre- 
quency at 200 Hz, and the current I D = ADI S is fed back 
to the signal circuit. The low-pass filter is used to phase- 
shift by 7r/2 the feedback current I D , which is then pro- 
portional to the derivative of the oscillating current of the 
mode. As a consequence, the effect of the feedback cur- 
rent is equivalent to that of a viscous damping, similarly 
to the case of a mechanical oscillator subjected to a force 
proportional to its velocity. This additional damping can 
be represented by an equivalent resistor R D in series with 
R k , which can be calculated from the model of Fig. [2] 
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In particular, in our experiment R D = \AD\uL in , as 
| AD | <C 1 and AD is almost purely imaginary. Induc- 
tance and capacitance of the mode are not significantly 
modified, provided that the feedback current phase-shift is 
close to 7r/2. Therefore, we write the total impedance of 
the circuit under feedback cooling conditions, for u ~ u k , 
as Z' k (to) = R k + R D -\-iujL k -\-l/ (iu>C k ). As customary, 
the relative strength of the feedback damping can be ex- 
pressed by the ratio of the feedback to the intrinsic damp- 
ing resistance (referred to the k-th mode) g k = R D /R k . 
Then, the quality factor of the k-th mode under feedback 
cooling is reduced to Q' k = Q k /(1 + g k ). 

According to the fluctuation-dissipation theorem, the 
thermal noise in the k-th normal mode is generated by a 
voltage noise source Vth with single- sided power spectral 
density Sy th = ^kBT R k , where T is the temperature of 
the thermal bath. This voltage noise is the effect of the in- 
teraction of the resonators with the microscopic degrees of 
freedom of the thermal bath, and therefore it is not affected 
by the feedback cooling. For to ~ uj k the power spectrum 
of the current noise induced in the resonator is: 
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The Lorentzian shape of the spectrum is determined by the 
feedback-reduced quality factor Q' k , while the prefactor is 
determined by the intrinsic quality factor Q k . Integration 
of Eq. ([3]l over frequency yields the total mean square cur- 
rent noise associated with mode k: 

k B T Q' k _ k B T 1 



I. 



(4) 



Lk Qk L k 1 + Qk 
so that, similarly to the quality factor, the effective temper- 
ature of the mode under feedback cooling is reduced to: 

To 



(5) 
1 + 9k 

As the temperature reduction is determined by the Q- 
reduction, the larger is the initial quality factor, the larger 
is the achievable feedback cooling. 

We can also rewrite Eq. (O in terms of the mode temper- 
ature instead of the bath temperature: 

4k B T k uJk 



\2 ' 



(6) 



Q' k L k (uj 2 - ul) + (u k u)/Q' k) 

This is precisely the expected power spectrum of a passive 
resonator with quality factor Q' k at thermal equilibrium at 
temperature T k . Actually, the resonator is not at thermal 
equilibrium at T k , but is rather in a dynamical equilib- 
rium between the thermal bath at T Q , and the measuring- 
feedback system, which acts as a very low temperature 
bath. 

In the above analysis we have neglected the back-action 
and additive noise of the SQUID amplifier, which are ex- 
pected to set a lower limit to the cooling efficiency. If the 
SQUID noise is taken into account, with the further sim- 
plified assumption of uncorrected noise sources, the fol- 
lowing refined expression of the mode temperature can be 
derived from the model in Fig. [2] 
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where S/ B is the spectral density of the SQUID measure- 
ment noise I n , fed back by the cooling loop, and SV„ is 
the spectral density of the SQUID back-action noise V n . 
For our present setup, the back-action contribution is al- 
most negligible with respect to thermal noise. According 
to Eq. (O, there is an optimum value of g k , for which the 
temperature achieves a minimum. This minimum achiev- 
able temperature is slightly dependent on the considered 
mode, and is of order 40 fiK for our present setup. 

We point out that feedback cooling of the modes does 
not improve the sensitivity of the system as GW detector. 
In fact, the cooling is due to a modification of the effective 
response of the system to any kind of excitation. Therefore, 
it suppresses in the same way both the thermal noise and 
the external signal originated by an impinging GW. 

To test the feedback cooling technique, we modified 
the standard operating conditions of the AURIGA detec- 
tor, and set four different values of the feedback cooling 
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FIG. 3: (Color online) Power spectrum of the current noise in 
the normal modes, as measured by the SQUID amplifier. The 
noise spectra are related to four different feedback settings. Each 
noise spectrum is well-fitted by a proper combination of three 
Lorentzian curves. Each Lorentzian peak is labelled by the cor- 
responding effective temperature, measured in mK. 



gain D, corresponding to four values of the relative feed- 
back damping g k . For each setting, we measured the power 
spectral density of the current detected by the SQUID sen- 
sor, and averaged for a time of roughly one hour. The 
four spectra are shown in Fig. [3] The overall power spec- 
tral density for a given setting can be accurately fitted 
by a proper combination of three Lorentzian curves, one 
for each mode. In fact, differently from previous exper- 
iments focused on a single mode, we simultaneously cool 
all three modes. The three-mode fitting curve is also shown 
in Fig. [3] superimposed on the corresponding noise spec- 
trum. For uj ~ Uk, the fitting function can be approxi- 
mated by the single-mode expression Eq. (O, from which 
we can infer the effective temperature of each mode T k . 
The experimental values of T k as function of the damping 
ratio (1 + gk) 1 are shown in Fig. [4] For a given feedback 
setting, different values of g k are associated to the three 
modes, because the intrinsic resistances R k are in gen- 
eral different. For instance, in our measurements g ± ranges 
from 190 to 2000 and g 3 ranges from 2200 to 30000. The 
values of gk are still small enough to make almost negli- 
gible the effect of the SQUID noise. Therefore, we expect 
the effective temperatures to follow the simple behaviour 
described by Eq. (0. The straight line in Fig.|4]represents 
T k for all 3 modes, calculated using Eq. (O without free 
parameters, with T fixed to the thermodynamic temper- 
ature of the bath, T =4.2 K. The data are in good agree- 
ment with the predictions of the model, even well below 1 
mK. The lowest achieved temperature are T\ = 2.0 mK, 
T 2 = 0.17 mK and T 3 = 0.20 mK. The lowest tempera- 
tures for modes 2 and 3 correspond to an average occupa- 
tion number (N k ) = k B T k /hwk ~ 4000. 

These results represent an improvement by more than 
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FIG. 4: (Color online) Effective temperature of the modes as 
function of the damping ratio. The straight line is the mode tem- 
perature predicted by Eq. lO, with To fixed to the value of the 
bath temperature T — 4.2 K, and no free parameters. The theo- 
retical limit at no feedback damping is also shown, to provide a 
graphical visualization of the achieved temperature reduction. 



the observation of any macroscopical quantum behaviour. 
In fact, the timescale r for decoherence in a resonator with 
quality factor Q and temperature T is proportional to to the 
ratio Q/T UM- 

Significant improvements with respect to the results pre- 
sented here are expected by performing a dedicated exper- 
iment on a mid-scale (0.01-1 kg) resonator at subkelvin 
bath temperature. Conventional dilution refrigerators can 
be used to cool kg-scale masses down to 10 mK ll2~lll . At 
the same time, resonators made of very low dissipation ma- 
terial, like silicon or sapphire, can reach quality factors as 
large as 10 9 lE2ll . Moreover, our capacitive-SQUID mea- 
surement system, similarly to SET-based ones [9], is nat- 
urally compatible with ultralow temperatures, due to the 
very low power dissipation, of order 10~ 10 W for typ- 
ical dc SQUIDs. Eventually, the maximum cooling ra- 
tio will be limited by the SQUID sensitivity, according to 
Eq. ((TV As state -of-art devices can approach the quantum 
limit 1123 , 24, 25|], resonator temperatures lower than 1 fiK 
are achievable at kHz frequency, corresponding to a single- 
digit occupation number of the quantum oscillator. 

We thank W. J. Weber for discussions and reading of the 
manuscript 



one order of magnitude with respect to the lowest tem- 
perature reported in literature for actively cooled macro- 
scopic mechanical resonators [6], and only experiments 
performed at much higher frequency [9] have achieved 
lower occupation number. Although in our experiment 
we actively cool electromechanical normal modes rather 
than purely mechanical modes, we believe that the com- 
parison is meaningful, as all three normal modes collect a 
significant fraction of the energy of the two mechanical res- 
onators. For instance, from the detector calibration we es- 
timate that the mode 2, the coolest one, collects about 36% 
of the energy deposited by an impulsive excitation in the 
bar resonator. Our cooling results are particularly relevant 
because of the enormous size of the elements of our system 
with respect to previously analysed micromechanical sys- 
tems. This overturns the common belief that the cooling 
should be easier for lighter resonators. To explain this ap- 
parent paradox, we observe that the the relevant parameters 
in determining the cooling capability are the displacement 
sensitivity and the intrinsic quality factor of the resonator, 
which sets the potential Q-reduction ratio. In our case, the 
latter is one order of magnitude higher than that usually 
achieved by micromechanical resonators. In fact, accord- 
ing to a well-established empirical rule, the quality factor 
of a mechanical resonator scales roughly with the volume 
to surface ratio, suggesting a limiting factor in the surface 
dissipation mechanisms [18]. As a consequence, large res- 
onators made of high Q material, with linear dimensions in 
the 1 cm - 1 m range, can easily reach quality factor as large 
as 10 6 -10 7 , whereas achieving Q > 10 5 in micron-sized 
resonators has as yet proven difficult II 1911 . We also notice 
that a very large quality factor is a necessary condition for 
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